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with similar onoe for the equatorial environment. In order to limit the 
analysis* this study will only consider the environment at 4GD km over the 
northern hemisphere during winter. Even with this limitation* the amount 
of information covered is still enormous. As o ronuLt* we hove further 
restricted tho study to perlnds of high color (sunspot number, R, of 1QQ) 
end Boomoflnetlc activity (flooinognotic activity level, Kp, of 6 0 ). Tho 
omghe^i o w-1-l l not be on the accuracy, but rather on the models nacosanry to 
adequately specify tho- shuttle environment. Motingo of tho actual module* 
data for other locations and conditions, and references to modulo not 
covered in tho report can bo obtained directly from tho author. 

Tho secondary object of tho study is to determine tho relative 
importance and sensitivity of different typos of onvi roamontol interactions 
os o function of tho environment. To accomplish this, whore possible* tho 
modollod environments hove boon used to predict tho level of the 
anticipated Interaction. Although this has proven to be a valuable output 
from the study* tho Interactions models ompluyod wore by necessity quite 
simplistic so that the absolute levels prodictod ere not Intended to bo 
accurate. Rather* tho results demonstrate potential parameter 
sensitivities and areas where the-environmentol models nosd to be improved. 

The report is organized into 4 sections dependent on tho environment 
being considered. In this study* only the natural atmosphere* geomagnetic 
field, ionosphere, and auroral environment at shuttle altitudes were 
considered. Models of the cosmic ray flux* radiation level, Solar 
electromagnetic flux* ambient electric field* gravity field, end debris 
environment will be presented at a later time. For each of the 
environments studied* an interaction Is modelled. For the neutral 

environment* the drag 16 computed. For the geomagnetic field, the Induced 
vxB electric field 16 estimated. For the ionosphere and auroraL 
environments, the vehicle to specs potential is estimated. The results of 
this analysis demonstrate* as would be anticipated* that there are indeed 
major differences in the environment between the equatorial and 
euroral/pola*’ environments that ere reflected in the interactions. 

THE NEUTRAL ATMOSPHERE 

By far the major env * ronmental factor at shuttle altitudes is the 
earth's ambient neutral atmosphere. Whether it be through drag or the 
recently discovered interactions with atomic oxygen, the effect of the 
r.eutrel atmosphere (predominately the neutral atomic oxygen) on spacecraft 
dynamics snd surfaces grsatly exceeds any of the other effects thet will be 
considered in this report. Currently there exi6t a number of models of the 
earth'6 neutral atmosphere. These models ere based on differing ratios cf 
data and theory. The 3 main sources of date at shuttle altitudes have been 
neutral mess spectrometers, accelerometers* and orbital dreg calculations. 
Without geiitg into detail, most models attempt to fit the observations with 
some fo.m of an algorithm thet includes the exponential fall off of the 
neutral density* tho effects of Increasing solar activity (particularly in 
the ultraviolet), the local time* and geomagnetic activity. Of these* the 
large variations associated with increasing gecmegnetic activity (end 
subsequent heating of the atmosphere) hove eluded adequate modelling by 
this fitting process. Unfortunately* it is clear from many sources (see* 
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for example, ref. 1) that those variations, particularly in density, over 
the auroral zone often dominate the neutroUonvIronmeni end that to date no 
adequate method of including these effects in. the models has been devised 
(some recent very sophisticated theoretical computer models do hold 
promise, however). 

With- the preceding caveat in mind, 2 models were used to compute the 
variations In drag due to the neutral- atmosphere at 400 km. These ere the 
Jacchia 1972 model (ref. 2) and the MSIS model (ref 6.- 3 end 4). These 
models are readily available in computer format and have been well 
developed over the last decade. For the purposes of this study, the 
Jacchia 1972 model results are presented (the MSlS model results deviate by 
about 20% from the Jacchia values on the average — a relatively small value 
given the much larger average uncertainties in the models themselves). 

Figure 1 illustrates the type of output that can be obtained, As stated 

earlier, the results are for the northern hemisphere (the reader is looking 

down on the north pole with the projection in. terms of equal latitude 

intervals) and 4Q0 km. The geomagnetic conditions for the model are for * 

FI 0.7= 2.2 x 1Q~ d W-nr-liz -1 (the 6olar radio flux at 10.7 cm; believed to 

be proportional to the extreme ultraviolet flux) and Kp=6 a . These give an 

exospheric temperature of about 1500 °K. ' 

Several features are apparent in the figure. First is the two-fold 
Increase in density from midnight to noon. Further, there is the 
pronounced shift by 2 hours of the peek in the density and temperature j 

maxima away from local noon. This well know phenomena results from the j 

rotation of the- earth and causes the peek in atmospheric hasting to occur i 

after local noon. The figure shows no clear feature associa-ted wi th the. 
auroral zone. This is directly due to the averaging techniques used in 
deriving models of this type which smooth out the density waves actually 
ob erved over the auroral zone. Even so, the model results ere useful in 
estimating the levels of atmospheric drag and, when the processes ace 
better knc^n, the Levels of shuttle "glow” and surface degradation. 

The major effects of the neutral atmosphere at 400 km result from the 
Impact of neutral perticle6 on spacecraft surfaces. Thi6 causes drag and 
surface damage. The standard expression for the drag force is: 

F(drag) = 1/2 pV 2 CD A = 

= —(300 - 5000) dynes (1) 

•IS q 

10 g/cnr at 400 km 
drag coefficient = 2.2 - 4.0 
cross-sectional area of spacecraft 
~50 m 2 (Frontal) for shuttle 
-400 nr (Bose) for shuttle 
spacecraft velocity 
7 .6 km/ o 

Comparing these values with Figure 1, it is evident that uncertainties in 
the orientation of the shuttle and lack of knowledge In the drag 
coefficient are equal to or greater than variations in the neutral 
environment at these Sltitudes. Given, however, the very real uncertainty 
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in the effects o-f auroral hosting, it is also apparent that these 
variations, i-f they ape greater than- a-f actor of 10 (which they can- ho), 

calculations 0 mo ^ or fton - tr<bu ^ or to uncertainties _1 n neu-tr-at drag 


MAGNETIC HELD- — 

Aside from the gravitational fiQld of the earthi the geomagnetic field 
at shuttle altitudes is the most accurately known. It can be crudelv 
modelled in terms of a tilted Ml 0 from geographic north) magnetic dipole 
of magnitude 8x10 G-cm . Numerous accurate, models of this field exist. 
Here we have used th . POGO model (refs. S end 6) as It is the basis of the 
International Reference Ionosphere (ZHl) model employed in the next 
section. This model is a straight forward expansion of fits to the earth's 
magnetic field in terms of spherical harmonics. The total magnetic field 
magnitude at 400 km according to this model is presented in Flg.s « 2. The 
surface field is seen to vary from a minimum of 0.25 6 near the ^uator to 
0.5 G over the polar caps. The existence of 2 peeks in the magnitude is 
reel and reflects the true complexity of the magnetic field in the 
auroral/polar cep regions (note: if vector components had been Included in 
. J J\? Ur ,° r iZ woUltl have baan obvious that the maximum at 270° east 
r. 9 U“tJ, 8 , th , 8 *" 9n,t,c Po'-e). Geomagnetic storm venations 

Mnnm VP J^f» l !? s ? hat , ‘ 91 9 60 that a,a ’' duMn3 o savers g.omegn.tic 
storm, magnetic fluctuations would be email compared to the average field — 

a marked contrast with the atmospheric end ionospheric environments! Even 
8 Brest complexity of the magnetic field over the poLes makes it 
difficult to use magnetic guidance systems in these regions— a fact long 
known to navigators. a 

Besides magnetic torques (which ere very system dependent) the 
£*£Sr* ,c f,ald ™ field in . l.rgTbo'yb'y the 


E = 0.1 (vxBJ V/m = 0.3 V/m 


( 2 ) 


where: 


v s spacecraft velocty = 
= 7.6 km/s 
B = 0.3 G 


lldoood b« a, t“h^. l °JLn^ U9 '’i , 15 " x 24 * » 33 «. potentlels of la V could bo 
Inducod f/olds .lUg^; £&££ ^ ° r Ur ° 0 

. , } n ^Bure 2, the Induced electric field for e vShicLe of ~.gno 

lli^V n 4 at *? n .i a8 been ealeulated * Aa would be anticipated, the largest 

ri:Hr^ 

,, „ " nor e , th 0 SB fields con reach values of nearly 100 mV/m (see ref 
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environment ae the particles have ambient energies of typicallv 0 1 aV (ram 

t ?° 1008 UkQ °* y0fin can roach severe! ev!' "hoover) aSf 21 
must be takeiMnto account -when -studying ionospheric fluxes. 

IONOSPHERE 

Given the importance of the ionosphere to radio and radar propagation. 
iftnlt n h ttrflr 8 i" 9 t0 f<wt th8t relatively, few models ere available for the 

L08 h SU ! f l P / i&1na 18 the fact tha * most of these models only 
#b h I densities— the most readily measureable quantity by ground 

® 8e08 8nd . . t * ,e roo6t important to radio propagation. The princioie 
Ionospheric model available based on observations is the International 

(r f* 8J * Th1e 18 th8 onl V readily wnmFiZSSS 
model that ; gives the electron and ion composition and temperature aa a 

function of longitude, latitude, altitude [65 to 1000 km), solar ecttvitv 

the mode! °i« oh® f U T 0t . r mber, r RJ| Md t1me { V B& r and local). Although 
the model is obviously limited (it is confined to R values of 100 or less 

whereas R values of 200 or greeter may occur during solar maxi mui) !! 

non.cheLa.fi 1. th. W .v.H.bl. ofimpV.h.n.1,. ^d^”f thfi1ol«XV.. 

In Figures 3 and 4, for the northern hemisphere, are presented severaL 

number^anei ou ^ put fro,n the IRI model. Figure 3 presents the electron 
number density end temperature at 400 km for R=100 in December. Unlike the 

S!iio a r^ am fJ r8tUre ’ th ® electron temperature Increases by a factor of 2^n 
illl 9 2 ! ! equator to the pole. Like the neutral density, however, the 

noon. tH ala0trOa d8nSUy 80810 18 ®h if ted by about 2 hours from local 

1ono8 P hepe * Primarily because of the corresponding high 
^ 8 * e !‘ neu t» , el oxygen. Is dominated by oxygen ions (4538 near Local 
midnight and below 30 latitude to 97* over the pole). Values for oxvaen 
ere presented in Figure 4. The temperature is assumed to be the sam^ for 
all ion species in this model (i.e., for 0 + , H*, He\ 0 P \ and N0 + ) and can 
not for physical reasons ever exceed the electron temperature 
Unfortunately, at 400 km for R=100 or larger, the IRl3ei 
occasionally predict Ion temperatures far in excess of the electron 

set^of^ata* tR<inoi pef ^ act8 th ® fact that the model is based on a. limited 
set ot data (R<100) and needs Improvement. Theoretical models axiet that 

oTiu'S t E\ l " theM " 0ltete — «<" «■ cunbereoBfi 1 to^^ueabta 
on all but the largest computers. 

ion Sul**? 5 » “ 8ln 8 8 8lm P Le 1 “dimensional, "thin sheath" ram model far 
ion collection (described in ref. 9), potentials for the case of no 

FigureVa n ° P hoto °t^tron current were calculated based on 

Figures 3 and 4. The spacecraft to space potential varied from -0 2 V at 
the equator to -0.7 V at the pole-in rough agreement with oSsarvetions 

charging should n-Tbe°e" concern (^e: 

ar^eHir^ toTtooT " 1tH aXP08Bd M0h P0tBntl8L B - f8C88 
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AURORAL ENVIRONMENT 


The most dramatic changes in the earth's environment at shtctie 
altitudes ere brought about by geomagnetic substorms.- Those changed- are 
reflected In visible auroral displays and In Intonse-particle and fisLd 
variations In the auroral region at shuttle altitudes. In thl6 section, a 

simple auroral flux model based on data provided by the Air Force 

Geophy6lc6_ Laboratory (courtesy M.. Smlddy and D. Hardys see papers by 
Smiddy end Hardy, this volume) 16 presented in order to estimate these 
effects. The data were provided- In the form of 7 sets of color contour 
plots of the electron number flux end energy flux in intervals of Kp from 0 
to 6. The plots were crudely approximated by a simple analytic function in 
geomagnetic local time and latitude and the geomagnetic Kp index. 
Although, the AFGL data were for about 800 km, no attempt has been made to 
correct for altitude In thl6 model. 

The crude model developed from the AFGL data was U6ed to. estimate the 
aurorel/poler cap electron temperature and number densities. These results 
for the northern winter hemisphere and a Kp of 6 Q are shown in Figure 6. 
They imply that there is a peak in the density of the auroral electron flux 
of about 1000 cm* 3 in the noon Sector while the auroral electron 
temperature Is 1 keV in the post-midnight sector. Although the validity of 
this crude resuLt will need to be compared with the actual AFGL data when 
they become available, the range of values 6hould at least be indicative of 
the characteristics of the average auroral fluxes, (comparisons with other 
data sources bear this out). 

The results. In-Figure 6 can bs-usad in conjunct-ion with the IRI date 
at 400 km to estimate the expected variations in spacecraft potential la- 
the auroral zone and over the polar cap6 (note: the auroral Ion fluxes 
should not contribute significantly to the ambient ion current So that 
their exclusion should not seriously aLter the results). When this 
Calculation was carried out, there wa6 little or no change from the results 
in Figure 5. Thi6 is not surprising as it is generally believed that the 
average auroral flux levels seldom exceed the ambient ion and electron 
ionospheric fluxe6. 

In order to estimate what auroral flux levels are in fact necessary to 
bring about significant increases in the spacecraft potential in the 
aurora l/po la r cap regions, the electron density and temperature in Figure 6 
were Increased by varying factors. Changes of a factor of 10 In either the 
temperature or density had little effect on the potential. A factor of 10 
in both the electron density and temperature did, however, bring about a 
significant increase In the potential — raising it from a few tehth6 of a 
volt negative to several thousands of volts in the early afternoon sector. 
These results are Illustrated in Figure 7. Such S large increase in the 
auroral flux may seem unrealistic but a careful review of auroral dote does 
imply that ocCassiohally intense fluxes 10 to 100 times that of the average 
flux may Indeed occur over narrow regions In the auroral zone (see, for 
example, Burke, this conference). 
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In carrying out the potential analysis, It. was found- that the details 
of the assumed charging model greatly affected the results. Specif icalLy, 
If a 1-dlmeneionaU thin sheath model was assumed, the auroral potentials 
could reach -8000 V when the Ion return current was equated to the cold 
ambient ion current. If the ion return current was assumed to be the ram- 
current, as was done here, the potential was about -1200 V maximum (-this la 
probably the more "realistic" assumption). If on the other hand- the ion 
return current- in the charging model was aasumsd to be for the- thiek 
sheath, orbit- 11 mi ted- case such os normally assumed at geosynchronous 
orbit, the potential was only -1 to -2 Vi This sensitivity to the details 
of the amount of return current is to be expected given the simplicity of 
the charging model and its resolution will need to await the development of 
more accurate charging models for the conditions at Shuttle altitudes. 


CONCLUSIONS 

This paper has brought together most of the elements needed to form a 
complete model of the ambient Shuttle environment for the purpose of 
studying spacecraft interactions. Emphasis has been on modelling, the 
interactions In the aurorel/polar cep regions where it was demonstrated 
that, although models of the average ambient environment (neutral 
particles, fields, Ionospheric particles, and aurorel/polar cap fluxes) ere 
probabLy satisfactory for many interaction study purposes, the intense 
variations in the auroral zone are not adequately modelled. These 
variations are known from in-situ observations to exlBt end to result in 
several orders of magnitude increaae In the charged particle fluxea and 
atmospheric heating which can similarly el ter the neutral compostion. It 
is only relatively recently that Long term statistical studies end-examples 
of extreme cases lave become available. It 1e to be anticipated that, in- 
the near future, models of the evironment will become increasingly 
sophisticated and capable of being used in modelling effects such as 
spacecraft charging much more occur tely than presented here. Even so, the 
results presented should assist current interaction studies in better 
assessing average levels of effects in the auroral/polar regions and in 
comparing equatorial and auroral/polar environments. The process of 
presenting the models has al6o clearly Indicated where improvements need-to 
be made in the existing models. Thl6 16 particularly true in the case of 
the auroral model due to the varying sensitivity of the principle 
Interaction to changes in the smbient environment li.e., spacecraft 
potential calculations). 

P. McConnell, M. Harel, and J. Slavin of JPL assisted in the 
collection and development of many of the models listed in this report. 
Any information on listings may. he obtained through them or the author 
directly. 


REEERENCES 

1. Straus, J. M.: Dynamics of the Thermosphere et High Latitudes. RfiVn. 
GbodHvb . . 16 1978, pp. 183-194. 

2. Jacchie, L. 8.: Atmospheric Models In the Region from 110 and 2000 km. 
CIRA 1972. Akedemi e-Verlao. Berlin, 1972, pp. 227-338. 


183 


3* Hedln, A. E. f Jr.; So La lx,. J. E.-; Evans, J. V»; Rebar, C. A.; Newton, 

6. P.; Spencer, N. W»; Kayaer, D. C.; Alcsyde, 0.; Bauer P.; Cogger, L.; 
end McClure, J. P.: A Global. Thermospheric Model based oh Maee Spectrometer 
end Incoherent Scatter Date, MSI S I , N2 Density end Temperatur a , J, 
Beophvn . Res.. 88 2139-2147, 1977a. 

4. Usd-in, A. E, Jr.; Rebar., C- A.; Newton, G. P.; Spencer, N. W.; 
Srlnton, H. C.; Mayr, H. G.; and Potter, W. E.: A Global Thermospheric 
Model based on Mase Spectrometer and Incoherent Scatter Data, MSIS 2, 
Composition, J. Geophvs.Res. .82 2148-2156. 1977b. 

5. Knecht, D. J.: The Geomagnetic Field. A1 r Force Surveys In Geophysics. 
No. 246. AFCRL-72-0570. 26 September 1972. 

6. Cain, J. C.; and Lancet, R. A.: The Geomagnetic Survey by the Polar 
Orbiting Geophysical Observations 0go-2 and 0go-4, 1985-1967, GSFC ReoJ<- 
612-68-502 . Greenbelt, Maryland, 1968. 

7. Foster, J. C.: An Empirical Electric Field Model derived from 
Chatanika Radar Date. jJL, Geophys. Res.. 86 19JB, pp. 981*987. 

8. Rawer, K. (Chairman): International Reference Ionosphere - IRl 79 
(edited by J. V. Lincoln and R. 0. Conkright). World Data Center A, Rot. 
UAG-82*. N ovember 1981. 


9. Garrett, H. 8.: The Charging of Spacecraft Surfaces. Rev. Geophys.. 19 
1981, pp. 577-616. 










Figure 2. - Polar view as in figure 1 for P060 magnetic field model. 
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(a) Density of auroral electron flux at 400 km. 

oure 6 - Polar view as in figure 1 of auroral flux model adapted from AFGL 
observations . Conditions correspond to Kp=6 0 and for day 357.5. 
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Figure 7. - Polar view as in figure 1 of spacecraft-to-space potentials pre- 
dicted for combination of IRI and auroral models. Potentials were computed 
as in figure 5 except that the auroral density and temperature from figure 6 
have both been multiplied by 10. 
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